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TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 

ACCURACY 

DURATION 

READ AROUND 
DESCRIPTION 



ADDITION 



Complex Number Operations (DO I or SAIJQl) 
Closed 

0,1,2,3 

Addition - exact 

>-39 



Division + 2 



-38 
,-38 



Modulus _+ 2 
Division 11 1/2 msec. 
Modulus 10 msec . 



Multiplication + 2 

Addition 3 1/2 msec. 

Multiplication 6 msec. 

Negligible 

If modulus is desired a closed square root routine must "be 

used in conjunction with this routine. The address of its 

first word, j, must be in location k in the form j x 2~^° 

when the complex number routine is read in using the D.0.X 

The Complex Number Routine will perform the operations of 

the L orders, 66, and 7 J orders for complex numbers. It 

assumes that the real part (Rp) and imaginary part (IP) of 

any complex number are stored in consecutive locations in 

the memory and reference is made to the address of the 

real part only. Operands and results must all lie within 

o 
the unit circle, i.e. for any number involved (Rp) + 

2 
(IP) < 1. The negative of the IP of any resultant is 

left in R^ so that conjugates may he easily formed by 

merely storing in the IP location of the resultant. 

Entry to the routine is made as follows: 

i-1 LA a. LB K 



i 50 
i+1 26 



50 



n 



where a-., b , are the addresses of the RP of the operands, 
c 1 is the address for the RP of the result, i is the 
sequence number in the main routine, and n is the address 
of the first word in the Complex Number Routine. 
A,B = 0,1,. ..,6 so that the L orders are interpreted in the 
usual way for real and imaginary parts separately, e.g. 



-2- 



MODULUS 



MULTIPLICATION 



DIVISION 



i-1 L5 a_ 1(0 b 1 

i 50 c 50 i 

i+1 26 n 

forms N(a x ) - N^) = N^) and N(a ) - N(b ) = N(c 2 ). 

L7 is -written for both LA and LB using the same address 

a., in each and is interpreted to form 



r 2 ^ 2 

a = J a., + b. 



"1 



1 



and place it in c. and in c « 

Multiplication is accomplished by entering the routine 

with the sequence 

i-1 50 a x 7 J b 



50 c x 50 i 



i+1 26 



n 



The complex product with RP = a.b - a~b and 
a p b + a b is placed in c 1 and c . 

Division is entered similarly except using L5 a 66 b in 

? 2 
i-1. This puts RP = i 3 -^ ^ + a 2 ^ )/(b + b ) in c and 

IP = (a^ - a^J/O^ + b 2 ) in c 2 - 
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S5 F 








1 


46 22L 
Ii4 52L 




Plant e* 




2 


42 24L 
46 23L 




Plant link 
Plant e 




3 


LO 53L 
42 4L 








4 


4l 3F 










L5 F 


By 3* 


HR = (.. a .« b ) 




5 


46 35L 
46 32L 




Plant a. 
n 




6 


46 45L 
42 26L 




Plant ^ 




7 


42 31L 
42 36L 




tt 




8 


42 45L 
42 4"7L 




n 




9 


L4 52L 
40 20L 




NE^ = ( .< a 2 ,. b 2 ) 




10 


46 30L 
46 37L 




Plant a 




11 


46 40L 




B 




• /"^' 


46 45L 




:« 




12 


42 2SL 
42 30L 




Plant bp 




13 


42 35L 
42 40L 




11 




lif 


42 44l 
42 4% 




n 
n 




15 


00 2QF 
46 28L 




KR 1 = (*. h 2 ..) 
Plant b p 




,16 


46 4% 




IT 






T.n k|pt. 




NEL = l..b *J 
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PAGE 2 1 


IT 


46 2TL 
k6 kSL 




Plant Ta- 




18 


56 25L 
LO kjL 




Test for + or (X or +) 




19 


52'^TL 
26 20L 




Test for add or modulus 




20 


00 F 
00 F 


By 9 ! 


Add 




21 


\0 3F 
L5 F 








22 


40 (e 1 )F 

L5 5F 


By 1 


Store resultant 




25 


40 (c 2 >. 
LI 5F 


-3^2 


Set-up conjugate 




2k 


22 24L 




■\ 






22 ( )F 


By 2 


Link"'*--, 




25 


LO 26L 
56 kOL 




Test for X or * 




26 


T^ OF 
50 (b^F 


From 18 
By 6 






27 


TJ (\> 
kO 3F 


By IT 






28 


50 (b 2 )F 
TJ (b 2 )F 


By 15 
By 12 






29 


Lk 3F 
ho 5F 




a * b 




50 


50 (a 2 )F 
TJ (\)F 


By 10 
By 12 






51 


ho OF 
50 (b^F 


By T 






52 


TJ (a x )F 
26 35L 


By 5 






55 


L4 F 









LOCATION 


ORDER 
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34 


S5 F 
40 F 






35 


50 (a 1 )F 


By 5 






7J (b 2 )F 


By 13 




36 


40 2F 








50 (b^F 


By 7 




37 


7J (a 2 )F 
LO 2F 


By 10 




38 


66 3F 
S5 F 






39 


40 3F 
22 21L. 






40 


50 (a 2 )F 


By 11 


From 25 




7J ("b 2 )F 


By 13 .. 




41 


40 3F 
L5 F 






42 


LO 3F 








40 F 




a x b 


43 


50 (a 2 )F 


By 11 






TJ.(-b 1 )F 


By 8 




44 


40 3F 








50 (b 2 )F 


By 14 




4 5 


7J (a 1 )F 
L4 5F 


By 6 




46 


40 3F 








22 21L 






47 


L6 102 3F 








50 (b 1 )F 


By 8 


From 19 


48 


7J (b 1 )F 
40 F 


By 17 




4 9 


50 (b 2 )F. 


By 16 


|b| 




7J (b 2 )F 


By 14 




50 


L4 F 
50 ^m, 







B2 



LOCATION 
51 

52 
55 



ORDER 

26 3k 
26 22L 
00 IF 
00 IF 
00 2F 
00 2F 



NOTES 
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B2 



N(*0 = Add of square root routine 



